Robust Identification of Developmentally Active Endothelial Enhancers in Zebrafish Using FANS-Assisted ATAC-Seq

Identification of tissue-specific and developmentally active enhancers provides insights into mechanisms that control gene expression during embryogenesis. However, robust detection of these regulatory elements remains challenging, especially in vertebrate genomes. Here, we apply fluorescent-activated nuclei sorting (FANS) followed by Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) to identify developmentally active endothelial enhancers in the zebrafish genome. ATAC-seq of nuclei from Tg(fli1a:egfp)(y1) transgenic embryos revealed expected patterns of nucleosomal positioning at transcriptional start sites throughout the genome and association with active histone modifications. Comparison of ATAC-seq from GFP-positive and -negative nuclei identified more than 5,000 open elements specific to endothelial cells. These elements flanked genes functionally important for vascular development and that displayed endothelial-specific gene expression. Importantly, a majority of tested elements drove endothelial gene expression in zebrafish embryos. Thus, FANS-assisted ATAC-seq using transgenic zebrafish embryos provides a robust approach for genome-wide identification of active tissue-specific enhancer elements

Transcription factor Hlx controls a systematic switch from white to brown fat through Prdm16-mediated co-activation

Browning of subcutaneous white fat (iWAT) involves several reprograming events, but the underlying mechanisms are incompletely understood. Here we show that the transcription factor Hlx is selectively expressed in brown adipose tissue (BAT) and iWAT, and is translationally upregulated by beta3-adrenergic signaling-mediated suppression of the translational inhibitor 4E-BP1. Hlx interacts with and is co-activated by Prdm16 to control BAT-selective gene expression and mitochondrial biogenesis. Hlx heterozygous knockout mice have defects in brown-like adipocyte formation in iWAT, and develop glucose intolerance and high fat-induced hepatic steatosis. Conversely, transgenic expression of Hlx at a physiological level drives a full program of thermogenesis and converts iWAT to brown-like fat, which improves glucose homeostasis and prevents obesity and hepatic steatosis. The adipose remodeling phenotypes are recapitulated by fat-specific injection of Hlx knockdown and overexpression viruses, respectively. Our studies establish Hlx as a powerful regulator for systematic white adipose tissue browning and offer molecular insights into the underlying transcriptional mechanism.The transcriptional co-activator Prdm16 regulates browning of white adipose tissue (WAT). Here, the authors show that Prdm16 interacts with the transcription factor Hlx, which is stabilized in response to beta3-adrenergic signaling, to increase thermogenic gene expression and mitochondrial biogenesis in subcutaneous WAT

Jmjd3-Mediated H3K27me3 Dynamics Orchestrate Brown Fat Development and Regulate White Fat Plasticity

SummaryProgression from brown preadipocytes to adipocytes engages two transcriptional programs: the expression of adipogenic genes common to both brown fat (BAT) and white fat (WAT), and the expression of BAT-selective genes. However, the dynamics of chromatin states and epigenetic enzymes involved remain poorly understood. Here we show that BAT development is selectively marked and guided by repressive H3K27me3 and is executed by its demethylase Jmjd3. We find that a significant subset of BAT-selective genes, but not common fat genes or WAT-selective genes, are demarcated by H3K27me3 in both brown and white preadipocytes. Jmjd3-catalyzed removal of H3K27me3, in part through Rreb1-mediated recruitment, is required for expression of BAT-selective genes and for development of beige adipocytes both in vitro and in vivo. Moreover, gain- and loss-of-function Jmjd3 transgenic mice show age-dependent body weight reduction and cold intolerance, respectively. Together, we identify an epigenetic mechanism governing BAT fate determination and WAT plasticity

MELK Promotes Melanoma Growth by Stimulating the NF-kappaB Pathway

Melanoma accounts for more than 80% of skin cancer-related deaths, and current therapies provide only short-term benefit to patients. Here, we show in melanoma cells that maternal embryonic leucine zipper kinase (MELK) is transcriptionally upregulated by the MAPK pathway via transcription factor E2F1. MELK knockdown or pharmacological inhibition blocked melanoma growth and enhanced the effectiveness of BRAFV600E inhibitor against melanoma cells. To identify mediators of MELK function, we performed stable isotope labeling with amino acids in cell culture (SILAC) and identified 469 proteins that had downregulated phosphorylation after MELK inhibition. Of these proteins, 139 were previously reported as substrates of BRAF or MEK, demonstrating that MELK is an important downstream mediator of the MAPK pathway. Furthermore, we show that MELK promotes melanoma growth by activating NF-kappaB pathway activity via Sequestosome 1 (SQSTM1/p62). Altogether, these results underpin an important role for MELK in melanoma growth downstream of the MAPK pathway

dagLogo: An R/Bioconductor package for identifying and visualizing differential amino acid group usage in proteomics data

Sequence logos have been widely used as graphical representations of conserved nucleic acid and protein motifs. Due to the complexity of the amino acid (AA) alphabet, rich post-translational modification, and diverse subcellular localization of proteins, few versatile tools are available for effective identification and visualization of protein motifs. In addition, various reduced AA alphabets based on physicochemical, structural, or functional properties have been valuable in the study of protein alignment, folding, structure prediction, and evolution. However, there is lack of tools for applying reduced AA alphabets to the identification and visualization of statistically significant motifs. To fill this gap, we developed an R/Bioconductor package dagLogo, which has several advantages over existing tools. First, dagLogo allows various formats for input sets and provides comprehensive options to build optimal background models. It implements different reduced AA alphabets to group AAs of similar properties. Furthermore, dagLogo provides statistical and visual solutions for differential AA (or AA group) usage analysis of both large and small data sets. Case studies showed that dagLogo can better identify and visualize conserved protein sequence patterns from different types of inputs and can potentially reveal the biological patterns that could be missed by other logo generators

Two Contrasting Classes of Nucleolus-Associated Domains in Mouse Fibroblast Heterochromatin [preprint]

In interphase eukaryotic cells, almost all heterochromatin is located adjacent to the nucleolus or to the nuclear lamina, thus defining Nucleolus Associated Domains (NADs) and Lamina Associated Domains (LADs), respectively. Here, we determined the first genome-scale map of murine NADs in mouse embryonic fibroblasts (MEFs) via deep sequencing of chromatin associated with purified nucleoli. We developed a Bioconductor package called NADfinder and demonstrated that it identifies NADs more accurately than other peak-calling tools, due to its critical feature of chromosome-level local baseline correction. We detected two distinct classes of NADs. Type I NADs associate frequently with both the nucleolar periphery and with the nuclear lamina, and generally display characteristics of constitutive heterochromatin, including late DNA replication, enrichment of H3K9me3 and little gene expression. In contrast, Type II NADs associate with nucleoli but do not overlap with LADs. Type II NADs tend to replicate earlier, display greater gene expression, and are more often enriched in H3K27me3 than Type I NADs. The nucleolar associations of both classes of NADs were confirmed via DNA-FISH, which also detected Type I but not Type II probes enriched at the nuclear lamina. Interestingly, Type II NADs are enriched in distinct gene classes, notably factors important for differentiation and development. In keeping with this, we observed that a Type II NAD is developmentally regulated, present in MEFs but not in undifferentiated embryonic stem (ES) cells

ATACseqQC: a Bioconductor package for post-alignment quality assessment of ATAC-seq data

BACKGROUND: ATAC-seq (Assays for Transposase-Accessible Chromatin using sequencing) is a recently developed technique for genome-wide analysis of chromatin accessibility. Compared to earlier methods for assaying chromatin accessibility, ATAC-seq is faster and easier to perform, does not require cross-linking, has higher signal to noise ratio, and can be performed on small cell numbers. However, to ensure a successful ATAC-seq experiment, step-by-step quality assurance processes, including both wet lab quality control and in silico quality assessment, are essential. While several tools have been developed or adopted for assessing read quality, identifying nucleosome occupancy and accessible regions from ATAC-seq data, none of the tools provide a comprehensive set of functionalities for preprocessing and quality assessment of aligned ATAC-seq datasets. RESULTS: We have developed a Bioconductor package, ATACseqQC, for easily generating various diagnostic plots to help researchers quickly assess the quality of their ATAC-seq data. In addition, this package contains functions to preprocess aligned ATAC-seq data for subsequent peak calling. Here we demonstrate the utilities of our package using 25 publicly available ATAC-seq datasets from four studies. We also provide guidelines on what the diagnostic plots should look like for an ideal ATAC-seq dataset. CONCLUSIONS: This software package has been used successfully for preprocessing and assessing several in-house and public ATAC-seq datasets. Diagnostic plots generated by this package will facilitate the quality assessment of ATAC-seq data, and help researchers to evaluate their own ATAC-seq experiments as well as select high-quality ATAC-seq datasets from public repositories such as GEO to avoid generating hypotheses or drawing conclusions from low-quality ATAC-seq experiments. The software, source code, and documentation are freely available as a Bioconductor package at https://bioconductor.org/packages/release/bioc/html/ATACseqQC.html

The coordinate actions of calcineurin and Hog1 mediate the stress response through multiple nodes of the cell cycle network

Upon exposure to environmental stressors, cells transiently arrest the cell cycle while they adapt and restore homeostasis. A challenge for all cells is to distinguish between stress signals and coordinate the appropriate adaptive response with cell cycle arrest. Here we investigate the role of the phosphatase calcineurin (CN) in the stress response and demonstrate that CN activates the Hog1/p38 pathway in both yeast and human cells. In yeast, the MAPK Hog1 is transiently activated in response to several well-studied osmostressors. We show that when a stressor simultaneously activates CN and Hog1, CN disrupts Hog1-stimulated negative feedback to prolong Hog1 activation and the period of cell cycle arrest. Regulation of Hog1 by CN also contributes to inactivation of multiple cell cycle-regulatory transcription factors (TFs) and the decreased expression of cell cycle-regulated genes. CN-dependent downregulation of G1/S genes is dependent upon Hog1 activation, whereas CN inactivates G2/M TFs through a combination of Hog1-dependent and -independent mechanisms. These findings demonstrate that CN and Hog1 act in a coordinated manner to inhibit multiple nodes of the cell cycle-regulatory network. Our results suggest that crosstalk between CN and stress-activated MAPKs helps cells tailor their adaptive responses to specific stressors